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Abstract 

We study theoretically the effect of ultraslow phase and group velocities in an anisotropic meta- 
material. The ultraslow phase propagation is induced by the hyperbolic dispersion relation. While 
the inherent physics underlying the slow group velocity are collective operations of the frequency 
and spatial dispersion. We show that a Gaussian wave packet exhibits simultaneous slow phase 
and group velocities which depend on the choice of incident angles and principal axis angles. The 
anisotrocpic metamaterial slab can be constructed and the ultraslow phase and group velocities 
can be measured experimentally. 
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I. INTRODUCTION 



The unusual properties of the electromagnetic waves in a medium with simultaneously 
permittivity e and the permeability were firstly introduced by Veselago about forty years 
ago ^|. The waves propagation in such media are quite counterintuitive. For example, the 
direction of energy flow for a plane wave is opposite to the direction of propagation. After 
the first experimental observation using a metamaterial composed of split ring resonators 
(SRR) P, Q|, the study of such materials has received increasing attention over the last 
few years. As noted earlier, both e and /i are necessarily frequency dispersive in LHM. 
Since the frequency dispersion is important, the superluminal 0, O, Q and the subluminal 
propagation p, Q, Q in the LHM takes place. 

It should be noted that we neglect bianisotropic effects despite the fact that they might 
be important for the characterization of some rings. While negative refraction is most 
easily visualized in an isotropic metamaterial, negative refraction can also be realized in 
anisotropic metamaterial, which does not necessarily require that all tensor elements of £ 
and /X have negative values . If we just consider the spatial dispersion, 

the superluminal group velocity can be expected in the AMM with hyperboloid dispersion 

AMM constructed by SRR is highly dispersive, both in spatial 



relation Il5|l. While the rea^ 



and frequency sense hence It is very necessary to extend the previous work 

and take the frequency dispersion into account. 

In the present letter, we will investigate the simultaneous slow phase and group velocities 
of wave packet in an anisotropic metamaterial. The subluminal phase propagation is induced 
by the hyperbolic dispersion relations associated with the AMM. We describe a modulated 
Gaussian wave packet incident on the anisotropic metamaterial, which demonstrates in 
a straightforward manner that the peak of the localized wave packet displays interesting 
smultaneous slow phase and group velocities. 



II. HYPERBOLIC DISPERSION RELATION 



It is currently well accepted that a better model is to consider anisotropic constitutive 
parameters, which can be diagonalized in the coordinate system coUinear with the principal 
axes of the metamaterial. If we take the principal axis as the z axis, the permittivity and 
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FIG. 1: The unit cell of the AMM is fabricated by printing rows of the SRRs onto large sheets of 
the dielectric substrate. The remaining volume is air-filled. For the sake of numerical calculations 
we choose the same parameters in the experiment carried out in Ref. jlJ ) 

permeability tensors have the following forms: 



e.{uj) 
ey{uj) 
e,{uj) 



(1) 





ti= ^ly{u) , (2) 



where Si and /Zj are the relative permittivity and permeability constants in the principal 
coordinate system {i = x,y, z). It should be noted that the real AMM constructed by SRR 
is highly dispersive, both in spatial sense and frequency sense So these 

relative values are functions of the angle frequency uj. A general study on the shape of the 
dispersion relation as function of the sign of these parameters has already been offered in 
Ref. |ll|. In this work, we are interested in the case of AMM with hyperbolic dispersion 
relation. 

Without loss of generality, we assume that the wave vector locate at the x — z plane 
{ky = Qy = 0). Maxwell's equations yield a scalar wave equation for E-polarized field . In 



+ + IQxQz = (4) 



free space, the accompanying dispersion relation has the famihar form: 

where kx and kz are the x and z components of the incident wave vector, u is the frequency, 
and c is the speed of hght in vacuum. We assume that there is an angle between the 
principal axis and the z axis. For the given polarization, the waves equation yield the 
dispersion relation in AMM as 

;^ 

C2 

where and represent the x and z components of refracted wave vector, a, (3 and 7 are 
given by 

a = {e^cos^ + EzSin^ if), 

(3 = {Ex sin^ + £z cos^ ip) , 

^x^zf^y 
1 

7 = sin 2(y9 — £3, sin 2^9). (5) 

^x^zf^y 

The material axes are normal to the surface of the rings. The cell structure is shown in 
Fig. H The metamaterial is formed by cutting the sheet into strips which are lined up 
with the appropriate periodicity. The permeability, /i^, can be approximated by the Lorentz 
model 

, ,2 _ 2 

f^z = l ^ ^ 77; , 6) 

where Umo, ^mp, and Vm denote the magnetic resonate frequency, plasma frequency, and 
damping frequency, respectively. The characteristic resonance frequency fmo = 10.08GHz, 
plasma frequency fmp = 10.56GHz and = O.lGHz (/ = u;/27r). Ignoring the metallic 
structure, Ey and fi^ assume the values of the background material which is dominantly air 
leading us to use Ey = fi^ = ^ [3- 

We assume here that the electric field is polarized along the y axis. We choose to fizi^^) < 
we find the corresponding requency contour is a hyperbola as shown in Fig. El where a plane 
electromagnetic wave is incident from free space into the AMM. We choose the z axis to be 
normal to the interface, the x axis in the plane of the figure, and the y axis out of the plane 
of the figure. Due to the angular dispersion, the transmitted wave components may refract 
at slightly different angles. The values of refracted wave vector can be found by using the 



4 



FIG. 2: The circle and the hyperbola represent the dispersion relations of free space and AMM, 
respectively. The frequency contour is rotated by an angle ip. The incident wave vector k is parallel 
to the group velocity vg in free space. Because of the anisotropy, the group velocity vg in the 
AMM is not necessarily parallel to the refracted wave vector q. 

boundary condition and hyperbolic dispersion relation Q|. The 2;-component of the wave 
vector can be found by the solution of Eq. which yields 



the choice of sign of ensures that light power propagates away from the surface to the +z 
direction. 

We now determine the angle of phase refraction. The incident angle of light in free space 
is 9 1 = ta.n~^[kx/kz] and the refraction angle of the transmitted wave vector in AMM can 
be found by 6x = tan'^lq^/qz]. From the boundary condition at the interface z = 0, the 
tangential components of the wave vectors must be continuous, i.e., q^ = k^- Thus the 
refracted angle of wave vector in the AMM can be easily obtained. When the interface is 
aligned an angle with the optical axes of the AMM, the hyperbolic dispersion relations will 
exhibit some interesting effects. 

To get a deeper insight into the slow phase velocity, we plot the frequency contour in 
Fig. 121 The magnitude of q varies as a function of its direction. When the refractive wave 
vector q is approximately parallel to the asymptotic line of hyperbola, q can be very large. 
The wave front travels in the AMM with the velocity of Vp = uj/q, so the ultra-slow phase 
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FIG. 3: The slow phase velocity in an AMM with different principal axial angles oi ip = (solid), 
tt/6 (dashed), 7r/4 (dash-doted) and vr/S (dash-dot-doted). The frequency of the incident wave is 
/ = 10.15GHz. 

velocity can be expected. The slow phase velocity in an AMM with different principal axial 
angles as shown in Fig. El 

It should be noted that a wave group can be formed from plane waves with different 
frequencies or form plane waves with different wave vectors. Thus the difference in the 
phase and group velocities can be caused by media that are dispersive or anisotropic 
The group velocity is a very important physical quantity because it identifies the speed of 
the maximum intensity of wave packet. The group velocity in anisotropic media can be 
defined as 



vg = —Re 



'dD/dq.; 


— Re 


'dD/dq; 




dD/dio 




_dD/duj_ 





(8) 



where D^q^, qz, uj) = aql + ^q'^ + '-yq^qz—^^^/c^, and are unit cartesian vectors. Because 
of the anisotropy and angular dispersion, the group velocity is not necessarily parallel to 
the wave vector q. Due to the resonance effect below the magnetic plasma frequency, the 
permeability functions undergo large changes with frequency, which results in dD/du ^ 
dD/dqx and dD/du ^ dD/dqz, hence the subluminal group velocity can be deduced. 

In Fig. m the group velocities in the AMM with different principal axes angle are plotted. 
It should be noted that the slow group velocity in the AMM is induced by the hyperbolic dis- 
persion relation. While the inherent physics underlying the slow group velocity are collective 
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FIG. 4: The subluminal group velocity in an AMM with different optical principal angles of 99 = 
(solid), 7r/6 (dashed), 7r/4 (dash-doted) and 7r/3 (dash-dot-doted). The frequency of the incident 
wave is / = 10.15GHz. 

operations of the frequency and spatial dispersion. 

III. SUBLUMINAL GROUP PROPAGATION 

To obtain a better physical picture of subluminal group velocity in AMM, a modulated 
Gaussian wave packet of finite width can be constructed. The field intensity distribution in 
free space is obtained by the Fourier integral and angular spectrum representation. Following 
the method outlined by Lu et al. let us consider a modulated Gaussian wavepacket is 
incident from free space 

/+00 
d?k^{\i - ko) exp[z(k ■ r - ic<;(k)t], (9) 
-00 

where c<j(k) = ck. we assume its Gaussian weight is 

f{k^^k;) = -eM-o'{kl + k% (10) 
vr 

where a is the spatial extent of the incident wave packet. We want the Gaussian wave packet 
to be aligned with the incident direction defined by the vector ko = k^ cosOie^ + k^ sin6'/e2, 
which makes an angle 9i with the surface normal. 
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FIG. 5: The temporal evolution of the Gaussian wave packet in the AMM. The center frequency 
is / = 10.05GHz with the incident angle Oj = — vr/G. The optical axis angle (f = it/3 and spatial 
extent of the incident wave packet is a = 10cm. The peak of localized wave packet travels from 
point 1 to point 2 during the propagation time At with group velocity vc ~ c/100. 

Due to the angular dispersion and the anisotropy, the transmitted wave components may 
refract at slightly different angles. When the Gaussian wave packet enters the AMM, it will 
no longer maintain Gaussian, but becomes a tilted wave packet. Matching the boundary 
conditions for each k component at z = gives the complex field in the form 

/ + 00 
d'k^fik - ko)T(k) exp[z(q ■ r - iuj{q)t], (11) 
-oo 

where uj{q) = cq and T(k) is the transmission coefficient. 

A close look at the tilted wave packet shows that the subluminal propagation of the peak 
is induced by the hyperbolic dispersion relation. Fig. shows a closes view of the field 
intensity distribution of the wave packet propagating in the AMM. In a propagation time 
At, the peak of the tilted Gaussian packet travels from point 1 to point 2. In Fig. El we set 
the center wave vector with a incident angles 6j = —it/ 6. We mark on each position of the 
peak of wave packet at each of the two times. In the At = 100ns, the peak of wave packet 
moves from (0, 0) to (2.45, 28.64). This propagating velocity corresponds to 28.75cm/ns, 
which is almost exactly the analytical group velocity vg ~ c/100 in Eq. (jHl), and the slow 
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group velocity is demonstrated theoretically. 



IV. DISCUSSION AND CONCLUSION 

In summary, we have discussed simultaneous slow phase and group velocities of light in 
an anisotropic Metamaterial. It should be noted that the slow group velocity has cornple tely 
different origin from those described in isotropic LHM or ultracold gas of atoms 18, 191 1. 
where the slow group velocity is caused by the frequency dispersion of the medium. In the 
case discussed here, the inherent physics underlying the slow group velocity are collective 
operations of the frequency and spatial dispersion. As far as we know, this kind of slow 
group velocity has not been recognized before. It should be mentioned, however, that the 
shape of Gaussian wave packet is distorted once it is incident into the AMM. It is shown 
that the simultaneously slow phase and group velocities can be measured experimentally. 

Acknowledgments 

This work was supported by projects of the National Natural Science Foundation of China 
(No. 10125521 and No. 10535010), and the 973 National Major State Basic Research and 
Development of China (No. G2000077400). 



[1] V. G. Veselago, Sov. Phys. Usp. 10 (1968) 509. 

[2] D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, Phys. Rev. Lett. 
84 (2000) 4184. 

[3] R. A. Shelby, D. R. Smith, and S. Schultz, Science 292 (2001) 77. 
[4] R. W. Ziolkowski, Phys. Rev. E 63 (2001) 046604. 
[5] J. F. Woodley and M. Mojahedi, Phys. Rev. E 70 (2004) 046603. 
[6] S. D. Gupta, R. Arun and G. S. Agarwal, Phys. Rev. B 69 (2004) 113104. 
[7] E. D. Gennaro, P. V. Parimi, W. T. Lu, S. Sridhar, J. S. Derov and B. Turchinetz, Phys. Rev. 
B 72, 033110 (2004). 

[8] G. DolHng, C. Enkrich, M. Wegener, C. M. Soukoulis, and S. Linden, Science 312, 892 (2006). 



9 



[9] I. V. Lindell, S. A. Tretyakov, K. I. Nikoskinen, and S. Ilvonen, Microw. Opt. Technol. Lett. 
31, 129 (2001). 

[10] C. G. Parazzoli, R. B. Greegor, K. Li, B. E. C. Koltenbah, and M. Tanielian, Phys. Rev. Lett. 

90, 107401 (2003). 
[11] D. R. Smith and D. Schurig, Phys. Rev. Lett. 90, 077405 (2003). 
[12] D. R. Smith, P. Kohnko, and D. Schurig, J. Opt. Soc. Am. B 21, 1032 (2004). 
[13] Z. M. Thomas, T. M. Grzegorczyk, B. I. Wu, X. Chen, and J. A. Kong, Opt. Express, 13, 

4737 (2005). 

[14] H. Luo, W. Hu, X. Yi, H. Liu, and J. Zhu, Opt. Commun. 254, 353 (2005). 
[15] H. Luo, W. Hu, W. Shu, F. Li, and Z. Ren, Europhysics Letters 74, 1081 (2006). 
[16] J. A. Kong, Electromagnetic Wave Theory (John Wiley Sons, New York) 1990. 
[17] W. T. Lu, J. B. SokolofF, and S. Sridhar, Phys. Rev. E 69 (2004) 026604. 
[18] S. E. Harris and L. V. Hau Phys. Rev. Lett. 82, 4611 (1999). 

[19] L. V. Hau, S. E. Harris, Z. Button, and C. H. Behroozi, Nature (London) 397, 594 (1999). 



10 



